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We present the results of three-dimensional simulations of quasar polarizations in the presence of 
pseudoscalar-photon mixing in the intergalactic medium. The intergalactic magnetic field is assumed 
to be uncorrelated in wave vector space but correlated in real space. Such a field may be obtained if 
its origin is primordial. Furthermore we assume that the quasars, located at cosmological distances, 
have negligible initial polarization. In the presence of pseudoscalar-photon mixing we show, through 
a direct comparison with observations, that this may explain the observed large scale alignments 
in quasar polarizations within the framework of big bang cosmology. We find that the simulation 
results give a reasonably good fit to the observed data. 

I. INTRODUCTION 

In a recent paper [l| it was shown that pseudoscalar-photon mixing in the presence of correlated intergalactic 
magnetic fields might provide an explanation for the observed large scale alignment of quasar polarizations at visible 
wavelengths [H-Q . The alignment is seen over cosmologically large distances of order Gpc 041] • In , the authors 
argue that it is very unlikely that the effect originates due to interstellar extinction. The effect is puzzling within 
the framework of big bang cosmology since we do not expect astrophysical objects to be correlated over such large 
distances. There is also direct evidence that other properties of these objects, such as, polarization position angles at 
radio frequencies or position angles of parsec-scale jets, do not show any large scale alignment @. Hence the effect 
docs not seem to be intrinsic to these sources and might arise due to propagation. Furthermore the propagation effect 
must not affect radio polarizations since these do not show any alignment. Pseudoscalar-photon mixing is a good 
candidate to explain the alignment of quasar optical polarizations since it is negligible at radio frequencies. Large 
scale correlations in the intergalactic magnetic field could lead to such an alignment over large distances [l| . 

In [l[ , it was assumed that the intergalactic magnetic field may be obtained by a simple cosmological evolution of the 
primordial magnetic field within the framework of big bang cosmology. It was assumed that visible radiation 

from distant quasars is unpolarized at the source. It acquires polarization due to its mixing with hypothetical low mass 
pseudoscalars in background magnetic fields (l5l - [26| . Such pseudoscalars arise in many extensions of the standard 
model of particle physics [2714351 ]. It was shown that the linear polarization angle of electromagnetic radiation from 
uasars separated over large distances becomes aligned due to correlations in the background magnetic field pj. In 
l| the authors restricted their attention to quasars that are aligned in one dimension, along the line of sight. In the 
current paper we present the results of a general three-dimensional numerical analysis of this problem, and show by 
a direct comparison with observations that such an effect can explain the observed large scale alignments in quasar 
polarizations. 

The mixing of photons with pseudoscalars has many interesting astrophysical and cosmological implications [H, |36| - 
5311 . It affects both the intensity and polarization of radiation. Furthermore it also changes its spectral characteristics 
[54j-l58[. Extensive laboratory and astrophysical searches of these particles have led to stringent limits on their 
masses and couplings |38l 1471. 159MT9H . Some astrophysical observations also appear to indicate a positive signature of 
pseudoscalar-photon mixing [53l . |80( 

The intergalactic magnetic field is not known very well. In most treatments of pseudoscalar-photon mixing through 
intergalactic space it is assumed that the background medium may be split into a large number of uncorrelated 
domains. The domain size is typically taken to be of order Mpc with the magnetic field in each domain of order 
nG [43l . l44l , l46j . The magnetic field in different domains points in random directions and is uncorrelated. Here we 
assume that the background magnetic field is of primordial origin (Tol - [T3 |. We also consider discretized real space 
consisting of a large number of cells or domains of equal size. Each domain is assumed to be cubical of side tq. The 
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j th component of the magnetic field, Bj(r), in real space is given by the discrete Fourier transform, 

where V is the volume in real space and bj(k) is the j th component of the magnetic field in wave vector space. The 
magnetic field is assumed to be uniform in each domain. We may write the two point correlations of bj(k) as, 

(bi(k)b*(q)) = 5 k , q P l0 {k)M{k) 

= 5 k , q a%{k), (2) 

where crfj(k) = Pij(k)M(k), k = \k\ and Pij(k) = — ^r^J is the projection operator. The function M(k) shows 
a power law behaviour, 

M(k) = Ak nB , (3) 

where rig is the power spectral index. The power law dependence of this correlation function would lead to correlations 
in the magnetic field in real space over large distances. This is in contrast to the assumption made in most treatments 
of pseudoscalar- photon mixing through intergalactic space [HI, 0, HH ■ The constant A in Eq. [3] is fixed by demanding 
that 

Y J <Bi{r)B l {r)>=Bl (4) 

i 

where we assume the value of about 1 nG [TH, HH, [H| for Bq . 

The value of the spectral index ns has been obtained by making a best fit to the matter and CMB power spectrum 
[8l| . This fit gives, ub = —2.37. A fully scale invariant power spectrum would correspond to the value ns = —3. 
We expect that the intergalactic medium may be turbulent on short distance scales. In this case we may expect the 
magnetic field to obey a Kolmogorov power spectrum for distance scales much smaller than the size of the system. 
A similar model is also applicable to the galactic medium, where both the plasma density and the magnetic field 
are known to show a Kolmogorov power spectrum (83l . I&ij . Hence it might be more appropriate to assume that rig 
has some dependence on k. Here we shall simply treat ns as a fixed parameter, independent of k. The spectrum 
may also have a low k cutoff, /c m ; n , which will imply a cutoff on correlations in real space for distances larger than 
r max = k~^ n . Here we shall assume that r max is larger than the size of our system, which we take to be of the order 
of a few Gpc. Hence we may simply set it equal to infinity. The large scale correlations induced among quasar 
polarizations crucially depend on this parameter. If this parameter is much smaller than 1 Gpc, then the mechanism 
proposed in [l[ may not work. Alternatively we may argue that the observation of large scale correlations in quasar 
polarizations might be an indication that r max > 1 Gpc. We emphasize that a magnetic field with such large distance 
correlations may be generated within the framework of the big bang model. The perturbations generated at the time 
of inflation have wavelengths extending up to the horizon in the current era. The Cosmic Microwave Background 
(CMB) anisotropics show significant power in the quadrupole, which implies correlations over the entire sky. The 
matter spectrum, however, gets significantly modified due to evolution after decoupling and correlations at such large 
distances get suppressed. However, the magnetic field need not evolve in the same manner as the non-relativistic 
matter. 

In our analysis we first ignore cosmological evolution and assume a flat space-time. We are interested in radiation 
from quasars located at redshifts of order unity. In this case our neglect of cosmological evolution would lead to an 
error of order unity, which is comparable to other errors, such as in the modelling of the intergalactic magnetic field. 
So far most of the literature on pscudoscalar-photon mixing in a background magnetic field has ignored cosmological 
evolution. We next also take cosmological evolution into account. 

A useful statistic to test for alignment of quasar polarizations is defined as follows [H, @ . We first define a measure 
of dispersion, d k , in the neighbourhood of the k th quasar, 

d k = — Vcos[2(V^ + Ai_> fc ) - 2^ fe )]. (5) 
n v f—' 

2 — 1 

Here ipi are the polarization angles of the nearest neighbours of the quasar at position k. We include a total of n v 
nearest neighbours. The nearest neighbours also include the polarization angle of the k th quasar. The concept of 
nearest neighbours requires a measure of distance. In a curved space we do not have a unique definition of distance. 
Here we use comoving distance as well as angular diameter distance for this purpose. When we include cosmological 
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evolution, wc avoid this problem by locating our simulated sources at the precise positions of the observed sources. 
The polarizations at two different angular positions are compared by making a parallel transport from their location 
to the position of the k th quasar along the great circle joining these points on the celestial sphere. This induces the 
factor Ai_>fe, included in the definition of dk- The parallel transport is required since we are comparing polarizations 
located at different points on the surface of the celestial sphere. If this parallel transport factor is not included then 
the resulting statistic is not invariant under coordinate transformations Jf|. We next maximize dk as a function of 
ipk ■ The resulting value of tpk is interpreted as the mean polarization angle at the position k and the corresponding 
maxima of the function in Eq. [5] gives an estimate of dk . The statistic may now be defined as 0, @ , 

1 n B 

SD = — Z, d *Lx' ( 6 ) 

s k=l 

where n s is the total number of sources in the data. A large value of Sd indicates a strong alignment between 
polarization vectors. In our analysis we shall compute this statistic theoretically and compare the corresponding 
values obtained by observations 0-0, @ in order get an estimate of the model parameters. 

A potential problem with our hypothesis that the alignment is caused by pseudoscalar-photon mi xing is discussed 
in Ref. [HI . It is observed that in most cases these quasars show negligible circular polarization [H, l86l |. So far this 
has been seen only in a small sample. Pseudoscalar-photon mixing predicts a larger circular polarization and hence 
may not consistently explain the data. The authors argue that if we assume the incident light to be natural white 
light, then due to decoherence the circular polarization is predicted to be zero, consistent with observations. However, 
circular polarization is obscrvationally found to be very small even with a broadband filter. In this case theory does 
indicate a significant circular polarization which is not in agreement with the data. Here we address this issue partially. 
We point out that there exist many other effects which may cause decoherence and hence reduce the predicted circular 
polarization. For example, in most treatments the medium is assumed to be uniform. Alternatively one assumes a 
large number of domains with different properties but each individual domain is assumed to be uniform. In reality the 
medium, at any reasonable length scale, shows fluctuations both in space and time. Since the background magnetic 
field is not uniform, the pseudoscalar particle produced would not have the same energy and momentum as that 
of the incident photon. This acts as another source of decoherence and will suppress the circular polarization. Wc 
show this explicitly by considering a model space varying magnetic field. We find that if the magnetic field shows 
rapid variations in space, then the circular polarization is greatly reduced in comparison to the linear polarization. 
Furthermore we argue that a time varying medium may also lead to suppression of circular polarization. 



II. PSEUDOSCALAR-PHOTON MIXING IN AN EXPANDING UNIVERSE 

In this section we give the basic equations of pseudoscalar-photon mixing in an expanding universe (87l - [89| . We as- 
sume the standard Friedmann-Robertson- Walker (FRW) metric, with the curvature parameter k = 0, using conformal 
time. The metric may be written as, 



a 2 
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Qiiv - I g q _ a 2 n | - a 2 V^, ( 7 ) 
V 

where a is the scale parameter and rj^ the Minkowski metric. 



A. Mixing with pseudoscalars 

The action of the electromagnetic field coupled to a pseudoscalar field <fi may be written as, 

S = f d^^—g [- - -^F^P* 



1 



" —2 ±2 



(8) 



Here g$ is the coupling of <f> with the electromagnetic field. We have also written an effective photon mass term. This 
corresponds to the effective mass ui 2 a~ 3 it acquires due to the medium, where lo p is the plasma frequency. In Eq. [8] 
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we have ignored any self couplings of the scalar field. The equation of motion for the scalar field may be written as, 

" ^ + V *X = ~^(a 2 E) ■ (a 2 B) + m 2 a 2 X , (9) 

where we have replaced <j> by In Eq. [5J B{ + Bi = ^eijka 2 F^ k is the magnetic field, Ei = a 2 F° l the usual electric 
field, B the background magnetic field, and B the magnetic field of the wave. The Maxwell's equations in curved 
space-time may be written as, 

V . {a 2 E ) = - 9 -±V X .{a 2 B) + p ] Ac, (10) 
a \ a I 

Vx(a 2 J E) + ^^=0, (11) 



with 



dr] 

d{a 2 B) 



dr] 



= 0, (12) 



V x (a 2 B) -^.%xx a 2 E) + 9 -± { a 2 B)^ + p ) A, (13) 
or] a a or] V a I 

V.(a 2 B)=0, (14) 

where A represents the vector potential. In the above equations we have assumed a uniform background magnetic 
field and approximated B + B « B as \B\ \B\. Here we also ignore the term containing derivatives of the scale 
factor, a, since the time scale for cosmological evolution is much larger in comparison to other times scales in the 
problem. Now we take the curl of Faraday's Law, and use the above equations to get the wave equation for E, 

- ^ + VV*) = ^(a-E) + (15) 
or] a a or]^ 

The time dependence of the resulting solution may be expressed as, 



(a E) oc cxp(— iuji]) = exp 



-VjJ 



dt 



(16) 



We also note that u(t)a(t) = luq, where ojq is the frequency at the current time. 

B. Mixing solution 

We choose the comoving coordinate system with the z-axis as the direction of propagation. In a particular domain 
the magnetic field is assumed to be uniform. Only the component of E parallel to the background magnetic field 

mixes with \- We also define A. = ^ ^ and replace (a 2 E) by ujA. We write the field equations of A\\ and \ as J 

(u 2 + d 2 )( A j)-M( A j)=0, (17) 



where, 

-^(a 2 B 



M =\ <M ^ 2> ■ ( 1S ) 
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Here, B± is the transverse component of B and u> is the frequency of radiation at the propagation domain. Due to 
expansion the observed energy today is redshiftcd. Hence we replace u> — > — . Finally the mixing matrix may be 
written as, 



We solve the mixed field equations in a manner similar to that in Ref . [2J, [24[ . We diagonalize the mixing matrix 
M by an orthogonal transformation, OMO T = Mp, where, 



where the angle 9 can be expressed as, 



0=( C °% «), (20) 

1 smO cos& 1 



tan2^^^. (21) 



We assume that the mass of the pseudoscalar (m^) is negligible compared to the plasma frequency (w p ). 



III. SIMULATIONS 



We first generate the magnetic field numerically in wave vector space. This is relatively straightforward since it 
is uncorrelated for different wave vectors. The projection operator implies that the component of the magnetic field 
parallel to k is zero. The two orthogonal components are uncorrelated. It is simplest to use polar coordinates (k, 9, </>) 
in wave vector space. Eq. [2] implies that for any wave vector fc, = and bg and b§ are uncorrelated. We may, 
therefore, generate these by assuming the Gaussian distribution, 

fQ>$(k),b^(k)) = N exp 

where N is the normalization factor. This represents an uncorrelated Gaussian distribution for the two components 
of the magnetic field in Fourier space, corresponding to the wave vector k. Once we have these two components we 
can obtain the three Cartesian components of the magnetic field in wave vector space. Next we use Eq. [T]to obtain 
the three Cartesian components of the magnetic field in real space. The Gaussian random variates are generated by 
using the Numerical Recipes (90j code gasdev. The discrete Fourier transform is obtained by using the Recipes code 
fourn. 

The optical polarizations are propagated, taking pscudoscalar-photon mixing into account by the procedure de- 
scribed in Ref. [79j . We need to propagate over a large number of magnetic domains from the quasar to the observer. 
We choose a fixed "external" coordinate system. The transverse component of the magnetic field in each domain 
is aligned at some angle to the fixed coordinate axes. In order to propagate through each domain we first rotate 
the coordinates so that the transverse magnetic field aligns along one of the coordinate axes. We then use standard 
expressions for pseudoscalar-photon mixing (79j in order to evaluate the correlation functions of the electromagnetic 
and pseudoscalar fields after propagation through the domain. We then rotate back to the fixed coordinate system. 
This procedure is repeated for propagation through each domain till the wave reaches the observer. 



b 2 e (k)+bl(k) 
2M(k) 



(22) 



IV. RESULTS 



In this section we present the results of our numerical analysis. We first consider the case of flat space-time, ignoring 
cosmological evolution. Next we give results for the case of expanding universe. We set the following parameter values: 
plasma density n e = 10 -8 cm -3 , frequency v = 10 6 GHz, coupling — 6.4 x 10~ n GcV^ 1 and mass = 0. We 
shall use the spectral index ns = —2.37 but will also explore a range of values of this parameter. We perform most 
of our simulations on a grid size of 1024 x 1024 x 1024, with the observer located at the center of the grid. The grid 
size is determined by the computing power available to us. 
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FIG. 1: The statistic So as a function of the distance between sources for different choices of the domain size tq (in Mpc) and 
the number of points on the grid No- Here we have set ub = —2.37. 
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FIG. 2: The number of nearest neighbours as a function of distance. The upper curve (dashed) corresponds to angular diameter 
distance assuming a vacuum dominated universe. The lower curve (solid) is obtained assuming comoving distance in a matter 
dominated universe. 

A. Flat space-time 

Here we present our results ignoring cosmological evolution. We assume 200 quasars distributed randomly in space. 
The resulting statistic Sd as a function of the distance among sources is shown in Fig. [1] Results are presented for 
four different values of the domain size tq — 4,8, 16,32 Mpc. The number of points on the grid are suitably scaled 
so as to keep the total volume fixed. In Fig. [T]the statistic is computed by including all the sources that lie within a 
particular distance from a source. This fixes the number of nearest neighbours to be included for any distance. The 
real data has a total of 355 sources. The relationship between the number of nearest neighbours and the distance 
for the real data is shown in Fig. [5] We find that for a distance of 25 Mpc, the number of nearest neighbours is 
approximately 1. It becomes close to 2 only at a distance of about 200 Mpc. In Fig. [3] we show the fluctuations in 
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FIG. 3: The fluctuations in the statistic Sd for different choices of the random seed. Here we have set tib — —2.37, ra = 8 
Mpc and the number of points on the grid Ng = 1024. 
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FIG. 4: Comparison of the simulation results with observed data. The dotted curve with large gaps shows the statistic obtained 
using observations. The remaining curves are obtained by simulations using parameters ns = —2.37 (dotted with small gaps), 
jib = —2.98 (solid), and a completely random magnetic field (dashed). Here we have set tg = 8 Mpc. 

the statistic Sd for a given choice of parameters ris = —2.37, tq = 8 Mpc and the number of points on the grid, 
Ng = 1024. Results are shown for four different seeds used to initialize the simulations. We see significant fluctuations 
and hence to compare with data it is better to take the mean over several different computations. Finally in Fig. 0] 
we compare the simulation results with observations. The simulation results are shown for a mean of five different 
initializing seeds choosing tib = —2.37 and —2.98. Here we have set tq = 8 Mpc and the number of points on the 
grid equal to 1024. We also show results for the case where the background magnetic field is taken to be completely 
random. We find that the results for correlated magnetic fields show reasonable agreement with data. However, those 
with a random magnetic field are systematically below observations. We also find a relatively mild dependence on 
the exponent ns- We conclude that pseudoscalar-photon mixing in a background magnetic field, described by Eqs. 
[T] and [21 gives reasonable agreement with observations. The small deviations from the data can be easily adjusted by 
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fine tuning the parameter values. 

B. Expanding Universe 

In this section we present our results for the case of an expanding universe. We consider 355 sources, to exactly 
match the number of sources observed in the real data. The positions of these sources are also taken to be the same 
as in the data. The medium parameters, i.e. the magnetic field strength, the plasma density are taken to be the same 
as in the earlier non-expanding universe. Here these values are assumed to be their comoving values. The observed 
frequency uj of the electromagnetic wave is assumed to be 2 eV. The comoving domain size is taken to be 15 Mpc. 
The Hubble constant is set equal to 2.133/i x 10~ 42 GeV, with h = 0.7. Furthermore we assume matter dominated 
universe. As we shall see our results for the expanding universe are comparable to those obtained for flat space-time 
and hence such details do not make a very large difference to our final results. 

We test our model against the data of the 355 sources observed by Hutsemekers et al. On placing the 

sources in our simulations at the same positions as in the data, we perform the propagation and calculate the Stokes 
parameters. We calculate the coordinate invariant statistics «Sx>, given in Eq. [SI for the observed polarizations and 
for the polarizations obtained from our simulations. The distribution of the linear polarization both for observed data 
and simulations is shown in Fig. [5] For the simulations we show results for rig = —2.37 and —2.95. 

In Fig. [HI we show the statistic Sd for the data as well as the simulations. The simulation results are shown for a 
mean of five computations. We see that for both jib = —2.37 and —2.95 the simulations show reasonable agreement 
with data. As mentioned above here we have set the domain size tq = 15 Mpc. It might be more reasonable to 
choose a smaller domain size. We find that as we reduce the domain size to tq = 12 Mpc, our statistic, Sd, becomes 
significantly larger than data. Hence we find even larger correlations in the optical polarizations. It may be possible 
to play with parameters to make our simulation results agree with observations in this case also. For example, this 
may be accomplished by decreasing the background magnetic field and/or the pseudoscalar-photon coupling constant 
<70. However we postpone a detailed fit to future work since it requires extensive numerical simulations. 

We may get some idea about the dependence of the statistic, Sd, on the domain size by performing simulations 
on a smaller grid. For this purpose we choose a grid size of 256 x 256 x 256. We again take all the 355 sources but 
place them at a smaller radial distance with their polarizations and angular positions taken to be same as that of 
real data. The resulting statistic, Sd, as a function of the number of nearest neighbours is shown in Fig. [7] for four 
different domain sizes. The results are shown after averaging over 40 different simulations in order to reduce the effect 
of fluctuations. We find that Sd shows an oscillatory behaviour as a function of the domain size. It increases as we 
increase the domain size from 10 Mpc to 14 Mpc. However with a further increase from 14 Mpc to 16 Mpc it starts 
to decrease. Similar trend is seen for a larger grid of 512 x 512 x 512. This suggests that for the range of parameters 
we consider, Sd shows an oscillatory dependence on the domain size. Furthermore the variation is significant but not 
very large. These results suggest that for the distance scale corresponding to the real data, the results may also show 
oscillations. Hence we expect to find a suitable fit to data even for smaller domain sizes with other parameters similar 
to what were chosen for the fit shown in Fig. [5] 

V. MIXING IN A MORE GENERAL BACKGROUND 

As mentioned in the introduction a potential problem with our hypothesis is that the observed circular polarization 
in quasars is found to be negligible 0, HfjJ. However our model predicts a large circular polarization. We argued 
in the introduction that our model of the background medium, i.e. the magnetic field and the plasma density, may 
be unrealistic. The medium is expected to show fluctuations in time and space at any length scale relevant to 
propagation over cosmological distances. Here we have assumed that the magnetic field and the plasma density are 
time independent. Furthermore we have assumed that the the magnetic field is uniform in a particular domain. In 
this section we determine how these assumptions affect our predictions for circular polarization. 

A. Space dependent magnetic field 

We first examine how small scale spatial fluctuations in the magnetic field affect final results. We will continue to 
assume a uniform plasma density throughout this section. As discussed above, the magnetic field is likely to have 
fluctuations even on the scale of a single domain. Here we assume a simple model of space varying magnetic field, 
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FIG. 5: The polarization histograms for the observed 355 sources (solid black curve) and for our simulations in an expanding 
universe. The red dashed curve and the green dash-dotted curve show results for ub = —2.95 and ns = —2.37 respectively. 



riven by, 



B = Bo (cos 2 a x + sin 2 a y) 



where, 



a 



2- 

T 



(23) 



(24) 



and A is the wavelength for the B field fluctuations. The rate at which the magnetic field changes with position z can 
be varied by changing the wavelength. Furthermore we point out that the field has been chosen such that its mean 
is not zero. The A and x field equations with this background can be written as, 



id z ) 









A 2 


-m i 


= o, 









(25) 



where, 



M 



2u 



2ui 



y - S± B cos 2 a ~^-Bo sin z a 



fB cos 2 a \ 
^B sin 2 a 

1u / 



(26) 



Here, we have approximated (oj 2 + d 2 ) ps 2uj{uj + id z ) [l^l and we have taken A as [A\ x + A2 y). These equations 
are same as those given in Section [II B| with the scale factor a = 1. 

Here we do not assume that the background magnetic field is changing sufficiently slowly so that we can choose 
basis vectors such that one of them points parallel to the transverse background magnetic field at any point z along 
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FIG. 6: The coordinate invariant statistic Sd as a function of the number of nearest neighbours for the 355 sources. The 
simulated results are shown by green triangles (tib = —2.37) and red squares (ub = —2.95), while the black crosses represent 
the observed data. 



the path of the wave. This approximation would be valid if A is very large but would break down for small A. Instead 
we directly solve these equations numerically in terms of the density matrix which is defined as, 

/ < At (z)A 1 (z) > <A\ {z)A 2 0) > <A\ (z)x(z) > \ 
p(z) = < Al{z)A 1 {z) > < A* 2 (z)A 2 {z) > < A* 2 (z) x (z) > \ . (27) 
\ < X *{z)A 1 {z) > < X *{z)A 2 (z) > < X*{z)x{z) > J 

The Stoke's parameters can be written as, 

I(z) =< Al^Aiiz) > + < A* 2 {z)A 2 {z) > (28) 

Q{z) =< AUz)A 1 (z) > - < A* 2 {z)A 2 {z) > (29) 

U{z) =< A\{z)A 2 {z) > + < A* 2 {z)Ai{z) > (30) 

V(z) = i(< A\{z)A 2 {z) > -A* 2 {z)Ai{z) >) (31) 

We assume that the wave is initially unpolarized and propagate in this medium over a distance of 10 Mpc. The 
resulting circular polarization |V|/7 and the linear polarization \/ Q 2 + U 2 /I is shown in Fig. [8] as a function of 
the wavelength A of the background magnetic field. The figure clearly shows that for small A, which corresponds to 
rapid fluctuations, the circular polarization is almost negligible as compared to the linear polarization. The value 
of linear polarization for small A is comparable to that at large A, where it starts to show large fluctuations. The 
circular polarization becomes comparable to linear polarization if the background varies slowly. Hence we clearly 
see that circular polarization is much smaller than linear polarization if the background magnetic field shows spatial 
fluctuations on sufficiently small length scales. 

Here we have presented results only for a single domain. We do not attempt the more ambitious task of integrating 
over all the domains. That would be much more computationally intensive in comparison to the calculation performed 
in this paper since it would involve solving a system of differential equations in each domain. Furthermore it would 
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FIG. 7: The coordinate invariant statistic Sd as a function of the number of nearest neighbours for 355 sources. The sources 
are located at the same angular positions as real sources but their radial distance is reduced so that we can use a smaller grid 
of size 256 x 256 x 256. Results are shown for four different domain sizes tg = 10 Mpc (red squares), 12 Mpc (green triangles), 
14 Mpc (blue inverted triangles) and 16 Mpc (purple circles), after averaging over 40 different simulations. Here we have set 
n B = -2.37. 



depend on the precise model we use for the magnetic field in each domain. In any case our analysis in this section 
clearly shows that there exist very reasonable model of magnetic field which would lead to highly suppressed circular 
polarization. 



B. Time dependent background magnetic field 

Another important issue that we consider is the possible time dependence of the background magnetic field. If the 
background magnetic field is time dependent, the pseudoscalar produced due to mixing with photons will not have 
the same frequency as that of the incident photon. Hence its correlation with photons may be significantly reduced. 
The circular polarization arises primarily due to reconversion of pseudoscalars back into photons. This generates a 
relative phase among the different components of the electromagnetic wave and hence leads to circular polarization. 
We see this explicitly by considering the evolution of the density matrix in a particular domain, given in Ref. [79| . 
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FIG. 8: The variation of linear (solid curve) and circular polarization (dashed curve) with respect to the wavelength A of the 
background magnetic field. 



reproduced here for convenience, 
<A n (z)AUz)> 



- < A||(0)Af|(0) > 1 + cos 2 29 + sin 2 29 cos[z(A - A A )} 



- < 0(0)0* (0) > sin 2 29 - sin 2 29 cos[z(A - A A )} 



< 



0(O)^*| (0) > sin 29 cos 29 - sin 29 cos 29 cos[z(A - A A )} 



<A ± (z)A* ± (z) > 
< A±(z)AUz) > 



i sin 29 sin[z(A^ — A A )] 



<Ax(p)A* ± (fl) > 
<A ± (0)A^(0)> 

< A ± (0)r(0) > 



cos 2 9e iFz + sin 2 9 e 1 Gz 



sin f cos 



(' 



A Fz 



A Gz\ 



(32) 
(33) 



(34) 



Here the vector A = E/w, ui is the frequency of the electromagnetic wave, is the pseudoscalar field and A\\ 
and A± respectively represent the components of A, parallel and perpendicular to the transverse component of the 
background magnetic field. These equations give the density matrix elements after propagating through one domain, 
assuming a space and time independent magnetic field and plasma density. The circular polarization is governed by 
the correlation < A±(z)A^(z) >. Let us assume that initially this correlation is zero. Then after propagating through 
one domain, it becomes non-zero only if initially the correlation < A±(0)cj>* (0) > is non-zero. We have assumed that 
at the source all the cross correlations are zero, i.e. the beam is unpolarizcd and the pseudoscalar field is zero. Hence 
circular polarization can be produced only if, after propagating through some domains, the correlation < A±(f>* > 
becomes non-zero. However here we argue that if the background magnetic field is time-dependent, the frequency of 
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the pseudoscalar field produced due to mixing would be different from that of the electromagnetic field. Hence its 
correlation would be significantly reduced. 

The pseudoscalar field equation may be written as, 

0_V 2 + m^ = 50 E-B. (35) 

where B represents the background magnetic field. Let us assume that the background magnetic field also undergoes 
fluctuations in time. Hence we may make a Fourier decomposition of this field. Here we assume that it contains only 
one significant component with frequency f2, i.e. \3\ ~ cos(fii). The frequency of the electric field is uj. It is clear 
from Eq. [35] that the frequency of the pseudoscalar produced has to be oj' = uj ± fi. 

Let us assume that the incident beam has a narrow spectral distribution, specified by the function A(uj), which 
denotes one of the components of the electromagnetic wave. We require correlations such as, < A(uj')<f>*(u}') >, where 
uj' = oj ± fl. Since </>(u/) oc A(uj)e~ lu t , we have 

< A(u)')<jf (a/) > cx < A{uj')A*{uj) > (36) 

It is clear that if uj' is sufficiently different from uj, this correlation would be vanishingly small. We expect this 
correlation to decrease rapidly with increase in the difference between the two frequencies. We point out that the 
background medium is expected to have fluctuations on all time scales. It is only after averaging out short time 
fluctuations that we expect it to show a somewhat smooth behaviour. Hence it may not be reasonable to assume 
that ft is close to zero. We point out that, although these temporal fluctuations may significantly affect the circular 
polarization, the mixing phenomenon would still generate significant linear polarization, as long as the mean value of 
the transverse component of the magnetic field is significantly different from zero. This is because the electromagnetic 
wave component parallel to the transverse magnetic field will decay into pseudoscalars and the perpendicular compo- 
nent will not. Hence we find that circular polarization may be significantly reduced compared to linear polarization 
if the background magnetic field shows sufficiently rapid fluctuations in time. 



VI. CONCLUSIONS 



We have shown by direct simulations, both in a static and expanding universe, that the hypothesis of pseudoscalar- 
photon mixing is able to explain the observed alignment of optical polarizations from distant quasars for standard 
parameter values of the intergalactic magnetic field and plasma density. The value of pseudoscalar-photon coupling is 
taken to be close to its observational limit. The simulation results show reasonable agreement with data both with the 
distribution of linear polarization as well as the coordinate invariant statistic. Given the uncertainties in the medium 
parameters we have not performed the exercise of finding the best fit to the data. In any case before attempting a 
detailed fit it is necessary to perform further calculations to conclusively establish that pseudoscalar-photon mixing is 
indeed responsible for the alignment. One potential problem with this explanation, as discussed in the introduction, 
is the relatively large circular polarization predicted by this phenomenon 85]. However, there are many sources 
of decoherence which need to be properly taken into account before making a firm conclusion on the strength of 
circular polarization. We have explicitly shown that if the magnetic field shows fluctuations over the distance scale of 
individual domains, the circular polarization predicted by this phenomenon is significantly reduced in comparison to 
linear polarization. We find this phenomenon even for a relatively large wavelength of fluctuations, of order 1 Mpc, 
of the background magnetic field. Furthermore the agreement we find with data is potentially encouraging and hence 
this proposal may be taken seriously. It is possible that even if this precise mechanism is found not to explain the 
data consistently, a suitable generalization might work. We postpone such issues for future research. 
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